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Abstract 

It is shown how pair production itself would almost certainly prevent the 
astrophysical formation of macroscopic dyadospheres, hypothetical regions, 
extending many electron Compton wavelengths in all directions, where the 
electric field exceeds the critical value for microscopically rapid Schwinger 
pair production. Pair production is a self-regulating process that would dis- 
charge a growing electric field, in the example of a hypothetical collapsing 
charged stellar core, before it reached 6% of the minimum dyadosphere value, 
keeping the pair production rate more than 26 orders of magnitude below the 
dyadosphere value, and keeping the efficiency below 2 x 10 _4 a/M/M©. 



1 Introduction 

Ruffini and his collaborators [Damour & Ruffini(1975), Ruffini( 1998a), Preparata, Ruffini, & Xue(l 
Ruffini et al.(1998), Ruffini & Xue(1998), Ruffini(1998b), Ruffini( 1998c), Ruffini et al.(1999a), 
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Ruffini et al.(1999b), Ruffini(1999a), Ruffini(1999b), Ruffini(1999c), Ruffini et al.(2000), 
Preparata, Ruffini, & Xue(2000), Ruffini(2000), Preparata, Ruffini, & Xue(2001), Bianco, Ruffini, 
Ruffini et al.(2001a), Ruffini et al.(2001b), Ruffini & Vitagliano(2002), Ruffini et al.(2002a), 
Ruffini et al.(2002b), Ruffini (2002a), Ruffini (2002b), Ruffini(2002c), Ruffini(2003), 
Chardonnet et al.(2003), Ruffini et al. (2003a), Ruffini et al. (2003b), Ruffini, Vitagliano, & Xue(20( 
Ruffini, Vitagliano, & Xue(2003b), Ruffini, Vitagliano, & Xue(2003c), Ruffini, Vitagliano, & Xue(: 
Ruffini et al. (2003c), Ruffini et al. (2004a), Ruffini et al. (2004b), Ruffini et al. (2004c), 
Ruffini et al.(2004d), Ruffini et al.(2004e), Bernardini et al.(2004), Bianco k Ruffini (2004), 
Corsi et al.(2004), Fraschetti et al.(2004), Bianco & Ruffini(2005a), Bianco & Ruffini (2005b), 
Ruffini et al. (2005a), Ruffini et al. (2005b), Ruffini et al. (2005c), Ruffini et al.(2005d), 
Bernardini et al.(2005), Xue et al.(2005), Ruffini et al. (2006a), Ruffini et al. (2006b), 
Vagenas(2006)] have proposed a model for explaining gamma ray bursts that pre- 
sumes the initial existence of what they call a dy ado sphere, a macroscopic re- 
gion of spacetime where the electric field exceeds the critical electric field E c = 
mlc 3 /(he) ~ 1.32 x 10 16 V/cm for microscopically rapid Schwinger pair produc- 
tion [Sauter(1931), Heisenberg & Euler(1936), Weisskopf(1936), Schwinger(1951), 
Nikishov(1970)]. The difficulty of producing these large electric fields is a prob- 
lem with this model that has not been adequately addressed. 

There are at least two strong reasons for doubting that such large electric fields 
can develop over astrophysical scales in all directions (i.e., over length scales much 
larger than the electron Compton wavelength, or much larger than the collision 
regions of individual charged particles). First, it would be very difficult to de- 
velop sufficient charge imbalance for macroscopic electric fields to produce signifi- 
cant numbers of pairs. Second, even if macroscopic pair production could somehow 
be achieved, I shall show in this paper that this process is sufficiently self-regulating 
that it prevents the electric field from achieving a value that would produce pairs 
at even 10 -26 that of dyadosphere models, assuming that the field is extended over 
a three-dimensional region at least as large as the Schwarzschild radius of a solar- 
mass black hole and that the field develops over a time that is at least as long as the 
corresponding timescale. I conclude that it is highly implausible that macroscopic 
dyadospheres can form in outer space, and therefore invoking them for models of 
gamma ray bursts is not at all likely to be viable. 

The first reason for being extremely doubtful of the existence of astrophysical 
dyadospheres is that it is very difficult for a large charge imbalance to develop over 
a macroscopic region astrophysically, because of the very high charge-to-mass ratio 
of elementary particles. For example, the ratio of the electrostatic repulsion to the 
gravitational attraction of two protons, each of charge q and mass m p , is the square 
of their charge-to-mass ratio in Planck units, which is 




(1) 



This implies that if one had a spherical object, such as a stellar core, with a positive 
charge-to- mass ratio Q/M greater than the inverse of the charge-to- mass ratio of 
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the proton, m p /q ps 9 x 10~ 19 , the electrostatic repulsive force on the protons at 
the surface would be greater than the gravitational attractive force, so such protons 
would most likely be ejected. (If the object had a negative charge, electrons of mass 
m = m e would be expelled if —Q/M > m e /q ps 4.9 x 10~ 22 , lower by the factor 
of the mass ratio of the proton and the electron, m p /m e pa 1836, so the maximum 
value of the charge of such an object would be even less if it were negative.) 

If one takes the mass-to- charge ratio of the proton to be a rough estimate of the 
maximum charge-to-mass ratio of an astrophysical object (or else gravity would not 
be strong enough to hold in the protons that make up the excess charge) , then using 
the formulas of the succeeding sections, one can readily calculate that at the surface 
of a spherical object of radius R and mass M, the ratio of the electric field value, 
E, to the critical field value of a dyadosphere, E c = m 2 c 3 /(hq) ps 1.32 x 10 16 V/cm, 



Therefore, if protons can be ejected from an astrophysical object whenever the 
electrostatic repulsion exceeds the gravitational attraction, then the electric field is 
constrained to be more than 13 orders of magnitude smaller than the critical value 
for a dyadosphere [if the mass is greater than 1.2 solar masses, which would be a 
conservative lower limit on any mass that could contract to 2GM/(c 2 R) ~ 1]. For 
a negatively charged object, the corresponding limit would be more than 16 orders 
of magnitude smaller than the critical value for a dyadosphere. 

Although it seems very unlikely to work, one might seek to evade the electrostatic 
expulsion of protons by postulating that they are bound by nuclear forces to an 
astrophysical object, such as a collapsing neutron star core. The critical electric 
field E c that Ruffini and his collaborators used to define the minimum value for a 
dyadosphere would give the electrostatic force on an electron or proton of magnitude 
F c = qE c = (m e c 2 ) 2 / (he) ps 0.00132 MeV/fm, whereas nuclear energies of the order 
of a few MeV over length scales of the order of a fermi would give nuclear forces of 
the order of a few MeV/fm, about 3 orders of magnitude larger than the electrostatic 
force of a minimal dyadosphere on a proton. Ruffini has proposed [Ruffini(2006)] 
that since the gravitational force is mainly provided by the nucleons of a collapsing 
neutron star core, with a tiny fraction of the nucleons in the form of protons that are 
bound by nuclear forces to the surface of the neutron star core, a sufficient fraction 
can give the electric field originating the dyadosphere. 

My suspicion is that such a configuration would be highly unstable to pieces of 
the charged surface breaking off and being ejected by the huge electrostatic forces 
on them. However, I do not have a firm result of when this would definitely happen, 
and therefore I was led to do the calculations below with the second mechanism 
(self-regulation) for preventing the occurrence of a dyadosphere. This, as we shall 
see below, will almost certainly discharge any growing electric field well before it 
reaches dyadosphere values. This occurs essentially because astrophysical length 
scales are much greater than the electron Compton wavelength, which is the scale 
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at which the pair production becomes significant at the critical electric field value 
for a dyadosphere. Therefore, the electric field will discharge astrophysically even 
when the pair production rate is very low on the scale of the electron Compton 
wavelength. 

In particular, I shall assume that a macroscopic astrophysical electric field cannot 
develop faster than the fastest timescale for bulk motions, which I shall assume to 
be the final gravitational collapse timescale for the smallest stellar mass that can 
collapse into a black hole, of the order of M ps 5 microseconds or greater. If 
a macroscopic electric field cannot develop much faster than this by astrophysical 
processes, the discharge from the pair production will then self-regulate the field to 
keep it from getting near dyadosphere values. 

The calculations below lead to the conclusion that it would be very difficult 
astrophysically to achieve, over a macroscopic region comparable to the size of a 
black hole or larger star, electric field values greater than a few percent of the 
minimum value for a dyadosphere, if that. The Schwinger pair production itself 
would then never exceed 10~ 26 times the minimum dyadosphere value. 

Since the idealized model below does give pair production at macroscopically 
significant rates (though more than 26 orders of magnitude below that of a dyado- 
sphere by the original definition given above), one might revise the definition of a 
dyadosphere to include any macroscopic electric field which gives macroscopically 
significant pair production. Then (assuming that sufficient charge separation can 
somehow be achieved by forces necessarily much stronger than gravitational forces, 
to evade the limitations discussed above), the arguments and model below would 
not exclude the possibility of such a revised concept of a dyadosphere. However, we 
shall see below that the much weaker amount of pair production in the model below 
does not seem to have nearly high enough energy efficiency for models of gamma 
ray bursts. Therefore in this paper, unless otherwise stated, I shall stick with the 
original definition of a dyadosphere, given above and in more detail in the following 
section. 

2 Pair production rate and definition of a dyado- 
sphere 

The Schwinger pair production [Sauter(1931), Heisenberg & Euler(1936), Weisskopf(1936), 
Schwinger(1951)] gives a rate M of electron-positron pairs per four-volume (per 
three-volume and per time) in a uniform electric field that is [Nikishov(1970)] 




(3) 
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where (using Planck units G = h = c = 47reo) m ~ 4.185 x 10 23 is the mass of the 
electron, q = y/a xs 0.08542 is the magnitude of the charge of the electron, and 

irm 2 i\E c , A . 

if one defines the critical electric field strength E c to be 
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E c = — ^ 1.323 x 10 ib V/cm. (5) 

Ruffini [Ruffini( 1998a), Preparata, Ruflini, & Xue(1998), Ruffini et al.(1998), Ruffini & Xue(19' 
Ruffini(1998b)] has defined a dyadosphere to be a region of spacetime in which the 
electric field E is greater than the critical electric field E c . Therefore, it is a region 
in which w < w c = ir, and a region in which the electron-positron pair production 
rate is 

4 — 7T 

M > M c = — w 10~ 93 « 4.7 x 10 48 cm" 3 s" 1 « 7.4 x 10 58 M" 4 . (6) 

47T 7T 2 © v ; 

The last number means that if one had a dyadosphere over a cube of edge length 
M & w 1.5 km and over a time M© ~ 5 microseconds, it would produce more than 
7 x 10 58 pairs, and the positrons would have a total charge ~ 1O 2O M . 

This shows that the minimum dyadosphere pair production rate, though only 
47re~ 7r in units of the electron Compton wavelength, is utterly enormous at macro- 
scopic astrophysical scales. This strongly suggests that dyadospheres will never form 
over macroscopic astrophysical scales. Indeed, the calculations below confirm this 
for a simple spherically symmetric model and show that under extremely conserva- 
tive assumptions, the pair production rate from astrophysical gravitational collapse 
is likely to be always less than 10~ 26 of that of a dyadosphere. (Under the plausi- 
ble but not rigorous arguments of the Introduction that the charge is ejected when 
the electrostatic repulsion exceeds the gravitational attraction, the upper limit in a 
region of the order of the size of a black hole would be trillions of orders of magni- 
tude weaker than a dyadosphere, a factor of more than lO 10 ", and hence completely 
negligible.) 



3 Crude upper limit of pair production rate 

One can make the following crude estimate of the upper limit of the pair production 
rate from any macroscopic astrophysical process that has an electric field E = 
nE c /w = nm 2 / (qw) extending along the field direction a macroscopic distance L ^> 
-kpair = 2m/ (qE) = 2w/{nm) (where L pair is the minimum distance along this field 
over which it is energetically possible to produce an electron-positron pair with 
total rest mass 2m), and lasting for a time scale T ~ L/c. We assume here and 
henceforth that since we are assuming a macroscopic electric field that extends over 
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a distance much greater than L pair , we may take it to be nearly homogeneous over 
the pair-production distance L pair , so that the pair production rate is given very 
accurately by Eq. (3) for a uniform field. 

This pair production over a time T and length L along the electric field direc- 
tion will lead to a charge density produced, per area perpendicular to the field, of 
magnitude a ~ qAfTL, which will discharge the field if it becomes comparable to 
E/(Air). Therefore, the average pair production rate is roughly limited by 

M=——^ E = m2 (7) 
An w 2 AnqTL Aq 2 wTL ' y ' 

Now if for the given T and L one defines 

v m 2 q 2 TL cTL 

X "^^" (8xlO- 1 o cm)2 » 1 ' ( § ) 

the inequality (7) implies that we w ~ X or w ~ In X — In In X . This then leads to 
the ratio of the actual pair production rate Af to the minimal dyadosphere rate Af c 
being 

^< ^ ( 10 " 8cm ) 2 «l ( 9 ) 

K X(hiX -InlnX) cTL V ; 

for any macroscopic time and length scales, which would give X ^> 7r 2 e 7r . 

This strongly suggests that no matter what the situation, if the time and length 
scales for the electric field obey 

7r 3/2 e 7r/2 

cT ~ L > l(r 8 cm. (10) 

mq 

then dyadospheres of such a macroscopic size will not form. 

In the following sections, we shall find a more precise upper bound on the pair 
production rate in an example with spherical symmetry. This example confirms 
the limit above and strengthens the evidence that pair production rates in macro- 
scopic processes must almost certainly be many, many orders of magnitude below 
dyadosphere rates. 

It must be admitted that the bound above is heuristic and is not a rigorous 
mathematical result. For example, if one just wants a hypothetical example of 
how a spatially infinite dyadosphere could develop from initial conditions without 
one, consider the following collision of two plane electromagnetic waves with the 
electromagnetic field tensor 

F = E [0(t - z)dx A d(t - z) + 0(t + z)dx A d(t + z)], (11) 

where 8 denotes the Heaviside step function. 

For t < 0, there are two incoming electromagnetic fields, with one front at 
z — t — —\t\ moving in the positive z-direction and the other front at z = — t = \t\ 
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moving in the negative z-direction. These fields are nonzero and constant behind 
their fronts, for \z\ > \t\, each with electric field E in the positive x-direction, 
but with the field at z < —\t\ that is moving in the positive z-direction having its 
magnetic field, also of magnitude E, pointing in the positive y-direction, and with the 
field at z > +\t\ that is moving in the negative z-direction having its magnetic field 
of magnitude E pointing in the negative y-direction. Because the scalar invariants 
of the electromagnetic field vanish for t < 0, there is no pair production then. 

However, for t > 0, the two incoming electromagnetic fields overlap for \z\ < t 
and have their magnetic fields cancel out there. In the absence of pair production, 
there would just be an electric field 2E in the x-direction in this region, which is 
infinitely extended in the x- and y-directions. Therefore, if 2E > E c , there would 
then be a dyadosphere that is infinitely large (in two directions) for t > 0, until the 
pair production reduced the field below the dyadosphere value. 

This hypothetical example circumvents the heuristic argument for the limit above 
first by not having the time scale T ~ L/c, but more importantly by having the 
incoming electromagnetic field (at t < 0) have Fourier components of arbitrarily 
short wavelength, in order to produce the sharp jump in the field from zero at 
\z\ < — t — \t\ to its nonzero values at z < —\t\ and at z > \t\. If instead the fields 
of the incoming plane waves rose sufficiently slowly (on the microscopic scale of the 
electron Compton wavelength, though this certainly allows what would be consid- 
ered a rapid variation on macroscopic astrophysical scales) to dyadosphere values of 
the electric field (though the equal magnetic field magnitudes in the perpendicular 
direction would prevent the field from actually being a pair-producing dyadosphere 
in each incoming plane wave when they are separate), then when the waves started 
overlapping to cancel the magnetic field and give pair production, the pair produc- 
tion itself would regulate the field from ever attaining dyadosphere values. There- 
fore, unless somehow it is possible for macroscopic astrophysical electromagnetic 
fields to develop significant components of microscopically large wavenumbers (and 
hence significant changes in the fields over length scales of the order of an electron 
Compton wavelength), it seems very unlikely that macroscopic dyadospheres will 
ever develop astrophysically. 

4 Schwinger discharge of a charged collapsing core 

In this section we shall analyze the pair production and discharge of an electric field 
produced by the collapse of a hypothetical charged sphere or stellar core, ignoring 
the processes discussed in the Introduction that would most probably cause almost 
all of the excess charge of the sphere to be electrostatically ejected. 

For a sphere collapsing in finite time, there is only a finite time for the electric 
field to be discharged by the Schwinger pair production process, so the discharge is 
never complete, but instead leaves a residual electric field at each moment of the 
collapse. We shall calculate an upper limit for this value and show that it is always 
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more than a factor of 18 less than the minimum value for a dyadosphere. Because 
the pair-production rate is exponentially damped by the inverse of the electric field 
value, the pair production rate itself never exceeds a value that is more than 26 orders 
of magnitude below that of a dyadosphere. (The factor of the order of 10 26 , which 
we shall derive below, comes mainly from the fine structure constant multiplied by 
the square of the ratio of the Schwarzschild radius of a solar mass black hole to the 
Compton wavelength of an electron, which is why this factor is so large.) 

The maximum electric field is smaller for cores that collapse into larger black 
holes (because the discharge time during infall is greater), so for a very conservative 
upper limit on the electric field, we shall assume that the black hole that forms has 
one solar mass. Of course, we expect that the minimum mass of a black hole that 
forms astrophysically has significantly more than one solar mass, so the correspond- 
ing maximum electric field would be weaker (by a logarithmic factor such that the 
maximum pair production rate would go essentially inversely with the square of the 
mass of the black hole). That is, we are assuming that one solar mass is a very 
conservative lower limit on the mass of any black hole that forms astrophysically in 
the present universe (as distinct from, say, primordial black holes that might have 
formed much smaller in the very early universe; our calculations will not apply to 
those, but they will have had by now a very long time to discharge and so would 
also not be expected to have significant charge today). 

The maximum electric field is also smaller the slower that the core collapses 
(giving more time for discharge), so again to get a very conservative upper limit on 
the electric field, we shall assume that the core falls in as fast as is astrophysically 
possible, which is free fall with zero binding energy. We assume that it is not 
astrophysically possible to have the spherical outer boundary of a collapsing core 
moving inward with a velocity so high that it would have come from an unbound 
configuration with nonzero inward velocity at radial infinity, or that at smaller radii 
it could have non-gravitational inward accelerations. 

Thirdly, the maximum electric field is of course smaller the smaller the initial 
charge on the collapsing core. To get the maximum electric field possible, we shall 
assume that the initial charge-to-mass ratio of the collapsing core is unity, the largest 
value possible for which the electrostatic repulsive forces do not overwhelm the 
total gravitational attractive forces on the entire core and prevent the core from 
collapsing (since we are excluding the possibility that the core is shot in from far 
away or is otherwise pushed inward by nongravitational forces, which is not at all 
astrophysically plausible). 

Now of course if the charge-to-mass ratio of the core really were unity, the core 
would not fall until there is some discharge. But merely to get a conservative upper 
limit on the electric field as the core collapses, we shall for simplicity assume that 
the initial charge equals the mass but ignore the electrostatic repulsion, so that the 
core nevertheless falls in at the rate it would from purely gravitational free fall from 
infinity, with no reduction of the free fall rate by the electrostatic repulsion of the 
charge. The gravitational effects of the electric field outside the core would also 
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reduce the free fall time, but we shall ignore this effect as well and simply take 
the external gravitational field to be given by the vacuum spherically symmetric 
Schwarzschild metric. 

First we shall make a rough Newtonian estimate of the discharge rate for a 
collapsing charged core. Next, we shall derive and give approximate solutions of 
the relativistic partial differential equations describing the process more precisely. 
Finally, we shall calculate the efficiency of the process for releasing the energy of 
the core. 

4.1 Approximate estimate of the discharge rate 

We consider a positively charged spherical core of mass M freely collapsing into a 
black hole, with the surface at a radius R(t). To avoid having the electric field dis- 
charged by plasma (a likely occurrence in most astrophysical situations), for the sake 
of argument we shall assume vacuum outside the core (except for the electromagnetic 
field and pairs produced by it, which we assume will not significantly modify the 
Schwarzschild geometry; including such modifications would reduce the electric field 
even further than the conservative upper limits we shall find below). Outside the 
surface of the core, pairs will be produced, with the positrons moving outward and 
the electrons moving inward. (Since the core is expected to be highly conducting, 
virtually all of its excess charge would be at or near the surface, so inside the core 
there would be a negligible macroscopic field and hence negligible pair production.) 
One can calculate [Page(2006a)] that the interactions between individual positrons 
and electrons produced outside the core are utterly negligible. E.g., the probability 
for each particle to annihilate with an antiparticle turns out to be less than 10~ 26 . 
Also, the number of particles produced per mode is very small, so we need not worry 
about particle degeneracy and Pauli blocking in the pair production region outside 
the core. 

As the electrons pass inward through the outer boundary of the core over time, 
they reduce the value of the charge of the core, Q(t), limiting the value of the electric 
field outside, E(t,r) ~ Q(t)/r 2 , under the assumption (which can be verified from 
the results) that at any one time the charge contributed by the pairs outside the 
core radius R(t) is a small fraction of Q(t), so long as one does not go to such a 
huge radius that it includes the outgoing positrons emitted over a large fraction of 
the previous infall of the core. 

We also make the assumptions, which will be verified later (at least for M <C 
1O 6 M ) that Q(t)«M once the core gets within a few orders of magnitude of the 
Schwarzschild radius 2M, and that the total energy that goes into the pairs is also 
much smaller than M (so that the core mass M stays very nearly constant). 

Because the pair production rate per four- volume, M = (m 4 /4n)w~ 2 e~ w , de- 
creases exponentially rapidly with w(t,r) = 7im 2 /qE m 7im 2 r 2 /qQ(t), and because 
the pairs produced decrease Q(t) and hence increase w(t,r) at fixed r, there will be 
a self-regulation of w(t,r) (described more precisely by the differential equations of 
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the following subsection, though we do not need this for the approximate estimate 
of this subsection). In particular, if we define 



z{t) 



w(t,R(t)) 



E(R(t)) 



ttE c 



7rm 2 R(t) 2 
qQ(t) 



(12) 



which gives the pair production rate Af(t) = (rrfi / A-n)z~ 2 e~ z at the surface of the 
core (where the rate is maximal at that time, since the electric field has an inverse-r 2 
falloff at greater radii r > R), the self- regulation will keep z(t) changing only slowly 
with t and R(t). Hence Q(t) will vary roughly in proportion to R{t) 2 . 

This means that the logarithmic rate of change of Q(t) will be roughly twice the 
logarithmic rate of change of R(t). If we do a Newtonian analysis of free fall from 

rest at infinity, we find that dR/dt = — ^2M/R, which makes the logarithmic rate 

of change of R(t) equal to — 1/(2M) times (2M/i?) 3 / 2 . As the core approaches the 
horizon at R = 2M, this Newtonian estimate of the logarithmic rate will approach 

— 1/(2M), so the logarithmic rate of change of Q{t) will approximately approach 

— 1/M. (Here we are ignoring relativistic corrections, but it turns out that they 
make only a very small difference.) 

Now we can calculate dQ/dt from the pair production rate (as a function of Q) 
and set it equal to —Q/M when R ~ 2M to solve for Q and hence the pair production 
rate. Because the logarithmic rate of change of Q is enormously less than that given 
above for a dyadosphere, the pair production rate must be suppressed by rather 
large values of w outside the core. Since w ~ nm 2 r 2 /qQ increases proportional 
to the square of the radius, the suppression will rather rapidly increase with the 
radius outside the core (which we are now taking to be at R ~ 2M). This means 
that the pair production rate will decrease roughly exponentially with the radius 
and will reach a value a factor of 1/e smaller than at the core surface itself (where 
w = z = 7im 2 R 2 /qQ ~ Anm 2 M 2 / qQ) roughly when w increases by unity. Since the 
logarithmic rate of increase of w with r at r = R is 2/R, the point at which the pair 
production rate will have dropped by the factor of 1/e will be at r—R pa Rj [2z] <C R, 
the last inequality coming because z ^> 1 from the fact that the pair production 
rate is much less than dyadosphere rates. 

With a roughly exponential decrease in the pair production rate with radius, 
at a logarithmic rate greater by a factor of 2z ^> 1 than the logarithmic rate at 
which the radius grows, the total pair production rate per time is roughly the pair 
production rate per four- volume at the surface of the core [J\f(R) = (m 4 /4:7r)z~ 2 e~ z ], 
multiplied by the effective volume where most of the pair production is occurring, 
which in this case is roughly the area AttR 2 of the core surface (or of the black hole 
formed by the collapsing core) multiplied by the radial distance r — R pa Rj (2z) 
out to where the pair production rate per four-volume has decreased by a factor of 
1/e. (For the pair production rate per coordinate time t, the relativistic correction 
that makes the proper radial distance greater than r — R is compensated in the 
approximately Schwarzschild metric outside the core by the relativistic correction 
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that makes the proper time smaller than At by what is precisely the same factor in 
the Schwarzschild metric.) 

When the number rate per time is multiplied by the charge — q of each ingoing 
electron, one gets at R ~ 2M that 

- « - q U{R)^- ~ —j—. (13) 

For the logarithmic rate of change of Q to be roughly — 1/M, we set dQ/dt equal to 
-Q/M = -nm 2 R 2 /(qMz) ~ -4vrm 2 M/(gz). This leads to the equation for z, 



z 2 e z ~ Z, (14) 



where 



and 



Z^.lgj .V«3.4X10V, (15) 
* S S « f^V « 3.396 x 10- (16) 

M , N 

(17) 

is the ratio of the mass M of the freely collapsing core to the mass M of the sun. 
One can solve Eq. (14) numerically for z when /i — 1 (M = M ) to get what I shall 
call z*: 

z* w 57.6. (18) 
Then for any |21n/i| <^ In A, one gets the approximate solution 

zw^ + -^-lnuw 57.6 + 1.93 In— . (19) 

^ + 2 M 

This means that the ratio of the electric field value E(2M) at the surface of the 
collapsing core (when it enters the black hole) to the critical electric field value E c 
for the definition of a dyadosphere is 

ggjfl = * w * _ ar ln„ Q Q546 _ Q 0M83 [n 

£ c z ^(^ + 2) M 

Thus this approximate estimate would indicate that the electric field outside a 
charged collapsing core is always less than about 5.5% of that of a dyadosphere, 
differing from that of a dyadosphere by a factor of at least 18. (We shall see later 
that the numerical solution of the ordinary differential equation governing the charge 
during the collapse, including relativistic effects, agrees with the result above to 
about three decimal places, so the result above is quite accurate.) 
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From this solution, one can see that the pair production rate at the surface of 
the core when it crosses the horizon at R = 2M is 



W R ) ~ ?W 2 *-* - 7^ = 7^ = ^ - 7 >< 10 ~ 27 ( ^-Y^ (21) 
v y 4tt * 47rA/! 2 4g 2 M 2 A/i 2 \ M 0/ 

which for M > M or /i = M/M Q > 1 is more than 26 orders of magnitude smaller 
than the minimum pair production rate M c of a dyadosphere. 

It is also of interest to compare the energy density pe of the electric field and 
the energy density p p of the particles produced with the mass-energy density p s of 
the stellar core when the core surface is crossing the event horizon. Of course, the 
density of the core will be affected by its infall and is likely to be inhomogeneous, 
but to get a quantity with which to compare pe and p p , let us say the core density 
is its total mass-energy M divided by the flat-space volume formula 4nR 3 /3 when 
R = 2M, which gives p s = 3/(327rM 2 ). Then since pe = E 2 /(8n), one may readily 
calculate that as the core surface crosses the event horizon, 

^^^!,1.4xlO-f^V- 933 . (22) 

The energy density of the particles produced peaks at r = 1.5R — 3M because 
of a trade-off between the growing electrostatic energy qQ(R~ l — r _1 ) gained by the 
positrons being accelerated outward by the Q/r 2 electric field and the 1/r 2 falloff in 
the number density. At this peak, one may readily calculate that 

Ps 81q 2 z 2 \M Q J V ; 

a factor of 16/27 times that of the electric field at r = R = 2M. Therefore, for the 
case Mc3x 1O 6 M in which significant discharge occurs, when the core collapses 
into a black hole, the electric field and pairs produced have a negligible energy 
density compared with the core itself, so they would not be expected to contribute 
significantly to the energetics of the core collapse or of any subsequent process like 
gamma ray bursts. 

Thus we see that the approximate algebraic estimate is that the pair production 
rate is never more than about 10~ 26 times the minimum amount for a dyadosphere, 
and the energy density of the electric field and of the particles produced is only a 
very tiny fraction of the collapsing stellar core energy density. In the next subsec- 
tion we shall indeed confirm from the solutions of the general relativistic differential 
equations that this algebraic estimate is indeed a good approximation for the maxi- 
mum electric field. This result implies that it is very unlikely that dyadospheres can 
form from the collapse of charged cores, even if somehow all discharge mechanisms 
are eliminated other than the pair production itself. 
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4.2 Differential equations for the discharge 



In this subsection we shall derive and give approximate solutions for the general 
relativistic differential equations for the discharge of the collapsing core. The differ- 
ential equations will be derived under the assumption that the tunneling distance for 
a pair to come into real existence, L pair = 2m/ (qE) = 2w/ (nm) ~ 10~ 9 cm, is much 
less than the astrophysical length scales for the collapsing core, which is a very good 
approximation. The differential equations will be solved under the approximation 
that this tunneling length is also significantly greater than the Compton wavelength 
of an electron, which implies that w = nm 2 / (qE) ^> 1 and hence that the pair 
production is mostly confined to a radial region r — R ~ R/w <C R that is much 
smaller than the radius R of the surface of the charged collapsing core. This latter 
approximation is less good but still leads to errors of only a few percent. 

After the charged particles are produced in pairs, with rms transverse momenta 



yqE/n = m/y/w, they will be accelerated by the radial electric field (electrons 
inward and positrons outward, under the assumption here that the core is positively 
charged). Each time a charged particle travels a distance m/(qE) parallel to the 
electric field, it will gain an additional kinetic energy equal to its rest mass. Thus 
it will very quickly accelerate to a huge gamma factor and move very nearly at the 
speed of light [with an asymptotic error at radial infinity for the outgoing positrons 
of l-v « w 2 /(2n 2 m 2 R 2 ) « 2.87x 10" 30 (M/M o )- L933 for positrons created when the 
core surface reaches the horizon at R = 2M] , and in very nearly the radial direction 
[with an asymptotic angular error of about \/w/(nmR) « 3.16 x 1O~ 16 (M/M )~°- 983 
when the core enters the black hole]. As a result, one will get effectively a null 
4-vector of positive current density (highly relativistic positrons) moving radially 
outward and a null 4-vector of negative current density (highly relativistic electrons) 
moving radially inward. 

It is most convenient to describe this current in terms of radial null coordinates, 
say U and V, so that the approximately Schwarzschild metric outside the collapsing 
core may be written as 



Now we can write the nearly-null outward number flux 4-vector of positrons 
as n + = n^<9v and the nearly-null inward number flux 4-vector of electrons as 
n = n^djj. Since each positron has charge q and each electron has charge —q, the 
total current density 4-vector is 



The radial electric field of magnitude E = Q/r 2 has the electromagnetic field 
tensor 




ds 2 = -e 2a dUdV + r 2 (U, V)(d6 2 + sin 2 6d(j) 2 ). 



(24) 



J = 



qn + — qn = qn + <9y — qn_<9u. 



(25) 



F = Ee 2a dU A dV 



Q 



e 2<J dU A dV, 



(26) 



2 



2r 2 
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where Q = Q{U, V) is the charge inside the sphere labeled by (U, V) (and which 
is a function only of these two null coordinates, because of the assumed spherical 
symmetry). 

Then from Maxwell's equations (essentially just Gauss's law here), we may de- 
duce that the null components of the current density vector are 

■v v -2e- 2 °Q tU +Q' V 

3 = qn + = 



4nr 2 Anr 2 
jU _ u _ +2e- 2CT Q,y _ z Qf_ { - 

3 ~ qn ~~ 4vrr2 ~ Airr 2 ' {Z '> 

Although of course the current density 4-vector field is conserved, the number 
flux 4- vectors n + and n_ of the positrons and electrons are not. Their 4-divergences 
are each equal to the pair production rate M when we can neglect annihilations, as we 
can here with the density of pairs being sufficiently small. When these 4-divergences 
are written in terms of the charge Q(U, V), one gets the following partial differential 
equation for the pair production and discharge process: 

Q, uv = -argrVW = exp j . (28) 

In covariant notation, with 2 OQ = — 4e _2o "(5,(yy being the covariant Laplacian 
in the 2-dimensional metric 2 ds 2 = —e 2cr dUdV and with □ being the covariant 
Laplacian in the full 4-dimensional metric, the differential equations for the pair 
production and discharge may be written as 

SvrgrW = ^ exp ( ^I^Illf] = 2 U Q = DQ - - Vr • VQ (29) 

or, explicitly in terms of the electric field E = Q/r 2 , as 

Mqr 2 N = 2<1 ^f exp = 2 a(r 2 E) = rU {rE) - r 3 E □ . (30) 

Since for the sake of argument in this section we assume that the positively 
charged particles at the surface of the collapsing core do not escape to the outside, 
and since there is no electric field inside to produce particles there, the boundary 
condition at the surface of the collapsing core is that there is no outward flux of 
positrons there, so that Qjj = at the core surface. The boundary condition at 
infinite radii is that we assume that there are no incoming electrons there, but only 
outgoing positrons pair-produced by the electric field at finite radius, so at radial 
infinity, Qy = 0. The boundary condition in the infinite past is that we assume 
that Q is as large as it can be and still have the core collapse gravitationally, so 
there we set Q = M but then make the idealized assumption of ignoring the effect 
of the electrostatic repulsion on the collapse of the core. 
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For details of the derivation of an approximate solution of the relativistic partial 
differential equations (28)-(30), see [Page(2006a)]. A brief summary of this solution 
is the following: 

Let 

be used as a time coordinate on the surface of the collapsing core with Schwarzschild 
(circumferential) radius r = R. For our assumption of the core surface freely falling, 
from an initial condition corresponding to being at rest at radial infinity, along radial 
geodesies in the assumed exterior Schwarzschild geometry, one may readily calculate 
that v is the negative rate of change of R with respect to the proper time r along 
the core surface, —dR/dr, and that it is also the magnitude of the three- velocity 
of the infalling core surface that an observer at rest at fixed r would observe when 
the core surface passes the observer. It runs from at radial infinity to 1 when the 
surface of the core crosses the event horizon, which is the limit of the validity of the 
approximate solution. 

The core surface has Schwarzschild time t, Schwarzschild radius r, tortoise radial 
coordinate r*, and proper time r given as explicit functions of v by 

2 

t = 2M[--v~ 3 -2V- 1 -ln{l-v) + ln{l + v)}, 
2M 



v 2 



r* = r 

AM 



+ 2M In - l) = 2M[v~ 2 - 2 In v + In (1 - v) + ln(l + v)}, 

(32) 



3v 3 



Next, choose null coordinates (U, V) in the exterior Schwarzschild geometry that 
are both equal to v on the collapsing core surface. These are defined implicitly as 
solutions of the following equations in terms of t and r* in the exterior region: 

t-n = AM[-\u- z -\u~ 2 -U- 1 +\nU -In (1-17)], 

t + r* = 4M[--V~ 3 + -V~ 2 - V' 1 -\nV + \n(l + V)]. (33) 

Then in the region exterior to the collapsing charged core (assumed to start with 
Q = Qo<MatR = oo, and yet with the charge having no effect on the geometry 
or infall rate, for the sake of argument to give a lower limit on the discharge), the 
approximate solution [Page(2006a)] gives 

Q(u, V) « ^f- {i - \ m [liVTTP + 1) - '-(ViTp - iu ^ ' 

(34) 
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Here 
where 

J = J(z(U)) =2z- 2V^e z erfc^ = l - 1 + - i^Ll + 0(z" 4 ), (37) 

and where z = z(U) obeys the ordinary differential equation 

= -8Sz/U 

d u 1 + U + ^/{i + uy -aus' 

At the core surface, U — V — v, one has Q(v,v) = nm 2 R 2 / (qz), so z = z(U) = 
z{v) = nm 2 R 2 /(qQ) = ixE c /E at the surface. For M < 3 x 1O 6 M , and for initial 
core charge Q = Qo large enough that significant discharge occurs, but not too large 
to prevent the core from collapsing, 

M Qo 

< T7 < !> ( 39 ) 



3 x 1O 6 M M 

a highly precise approximate solution of Eq. (38) for z(v) is [Page(2006a)] 

3 45 
z(v) « g(v)-2\ng(v)+g(v)- 1 (A\ng(v) - -) + g{v)~ 2 {A In 2 g(v) - ll\ng(v) + — ) 

+ g (v)- s (™hi s g(v) - 30\n 2 g(v) + f\ng(v) - ^), (40) 



where 

g(v) = f(v) + In 



v 5 z 1 (v) 2 e Z1 ^ 5(1 + v) 5(l + 2u)(3-2u) 
+ ZJ{ Zl {v)) + Az 2 {v) 16z 2 {v) 2 



(41) 



Ai\m 2 M 2 

^ s (42) 

»(„) 5 /M - 2 h / M + 4 h *»> - ^ + L25 % 4 '" 2 ^ - y - 2 - 5 "» . 

(43) 

M (44) 

A less precise but still reasonably good, and much simpler, approximate solution 
of the ordinary differential equation (38) is given implicitly as the solution of 

9 , fAnm 2 M 2 \ 2 fAirm 2 M 2 \ q 2 m 2 M 2 , A . 

z e " [i^r) exp l^^j + (45) 
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or, equivalently term for term for Q = Q(R) at the core surface, 



7T 



EM 2 



Q 



exp 



irE c R 2 

. Q 



ttE c R 2 



cxp 



'ttE cj R 2 \ q 2 m 2 R 2 (IT 
Q + Air V 2M' 



(46) 



When f = ^2M/R is sufficiently small (or i? is sufficiently large) that the first 
term on the right hand side greatly dominates, little discharge has occurred, and 
Q(R) ~ Qo- When v gets sufficiently large, the second term on the right hand side 
greatly dominates, and then the discharge and self-regulation of the charge becomes 
important. This will have occurred by the time that the core reaches the horizon 
(v = 1) if M < 3 x fO 6 M and if Q /M > M/(3 x 1O 6 M ) so that there was 
initially enough charge for the electric field to get high enough during the collapse 
to lead to significant discharge. Then taking z 2 e z « q 2 m 2 M 2 /ir at the horizon leads 
precisely to the Newtonian result of Eq. (14) and maximum pair production rate 
~ 7 x 10~ 27 (Mq/M) 2 times that of a minimal dyadosphere rate with E = E c . That 
is, the approximate solution of the general relativistic partial differential equations 
for the discharge leads very nearly to the same result as the approximate Newtonian 
estimate. 

Using the more precise approximate solution given by Eq. (40) , or the numerical 
solution of the ordinary differential equation (38) (both of which agreed to several 
decimal places), a more precise formula for the ratio of the maximum pair production 
rate at the horizon to that of a dyadosphere was given in [Page(2006a)] as 



N{M) 



6.612 x 1(T 27 (^f 
V M 



1 + 0.000555 In 



M_ 



0.000018 In 2 



M_ 



.(47) 



Therefore, assuming that M = M Q or /x = 1 is a very conservative lower limit on 
the mass of a core that can collapse into a black hole, we indeed see that even if one 
can somehow start with Q = M when the core is very large, and somehow not have 
the charge on the core itself directly ejected by the enormous electrostatic forces 
(other than the discharge by the pair production process), the pair production rate 
would always be more than 26 orders of magnitude smaller than that of a putative 
dyadosphere. In particular, the minimum dyadosphere value would be more than 
1.5 x 10 26 times larger than the actual value. 

Another quantity of interest is the maximum value of Q/M of the black hole 
when it initially forms. For M 10 6 M Q , this is then 



Q_ 

M 



Anm 2 M Aiim 2 M 



qz(l) 



1 - 



qz* 



21n// x 

z* + 2, 



4.09 x 10" 7 



M 



Mr 



1 - 0.0336 In 



M 







M 



.(48 



0, 



This implies that stellar mass black holes would always have a very low charge-to- 
mass ratio, even if the star somehow had a large charge before the collapse and 
ensuing discharge by pair production. 
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Furthermore, if one were to envisage a collapsing core much larger than stellar 
masses, sufficiently massive that it could collapse into a black hole without discharg- 
ing significantly, it is very hard to imagine how the positive charge (e.g., protons) 
could avoid being electrostatically ejected. Even if nuclear forces were somehow 
effective in accomplishing that Herculean feat for neutron star cores, it would seem 
even much more unlikely that one could form a neutron-star-like core of very many 
solar masses, so that nuclear forces on the protons could conceivably be effective in 
overcoming the huge electrostatic repulsion if there were a significant charge imbal- 
ance. 

Therefore, I would actually be surprised if any black holes of astrophysical masses 
ever form within our universe (or our pocket universe with our values of the masses 
and charges of the electron, proton, and neutron) with values of their charges at all 
near their masses. To put it more concretely, I would predict that no astrophysical 
black hole ever has a detectable change in its geometry given by the energy in its 
macroscopic electric field. In other words, Q/M would always be so far below unity 
that the metric of an astrophysical Reissner-Nordstrom or Kerr-Newman black hole 
would be indistinguishable from a Schwarzschild or Kerr black hole. 



5 Energy efficiency of the pair production 

We have calculated that even with idealized conditions of a collapsing stellar core ini- 
tially somehow having Qq = M and somehow keeping its excess protons from being 
driven off by the extremely large electrostatic forces (though admittedly smaller than 
the nuclear forces within a nucleus), one cannot get the electric field to become large 
enough to produce pairs at a rate per four-volume within 26 orders of magnitude of 
the minimal dyadosphere rate. However, we did get astrophysically significant pair 
production during this idealized process (which might count as a dyadosphere under 
a revised definition of the term), and so one might ask what the energy efficiency of 
this process is, what fraction e of the mass-energy M of the stellar core is converted 
into outgoing positrons, to give them total energy eM. Here we shall show that 
unless M > M , the efficiency is very small, e < 1.86 x 10" 4 ^M/M & < 1. 

In this approximately Schwarzschild metric (24), as written in terms of radial null 
coordinates U and V, if a positron is produced at some (U, V) and then is accelerated 
to very high gamma factors essentially along the outward null line U = const., then 
the kinetic energy (as measured at radial infinity) gained during this acceleration 
by the electric field Q(U, V)/r 2 (U, V) is 

mV) -Jn —-Jv 2r*(U,V>) ' (49) 

where the integral is taken along the outward radial null line U = const, along 
which the positron approximately travels from its creation point at (U, V) to radial 
infinity at (U, oo). Almost all of the positrons are created sufficiently deep in the 
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electric field that their final kinetic energies, at radial infinity, far exceed their rest 
mass energies, so I shall ignore the latter. Then one may multiply the positron 
production rate by the energy gained by each positron and integrate over all of the 
exterior region where the production is occurring to get the total energy emitted 
during the core collapse as 

eM = — f £{U, V)Q )UV dUdV. (50) 
q J 

(The negative sign comes from the decrease of Q during the discharge, so Q : uv < 0.) 

The crucial quantity for estimating the efficiency is the radius R t where the 
self- regulation of the charge starts becoming important, say at the point where 
both terms on the right hand side of Eq. (45) or (46) become equal. For general 
\i = M/Mq and £ — Qo/M, the critical radius where the self-regulation starts 
becoming significant is [Page(2006a)] 

Rt = 2M/v 2 t » 5300 ^0-5049^0.5082 km (51) 

Thus if one starts with a hypothetical charged core collapsing freely from a very 
large radius into a black hole with initial charge Qo ~ M, the self-regulation of the 
charge will start to become important when the the core gets to a radius of the 
order of the radius of the earth, and that is when there will start to be significant 
pairs being produced. The proper time left during the free-fall collapse after the 
core surface crosses this radius is then [Page(2006a)] 

At = -Mjv\ w 0.499 /i a25? £ a762 seconds. (52) 
3 

Therefore, if Qo ~ M, the burst of positrons that will be emitted as the core collapses 
from R = R t to R = 2M lasts a time that is of the order of one second. 

Then when one evaluates the integral in Eq. (50), one gets [Page(2006a)] that 
the efficiency of the conversion of the core mass M into outgoing positrons of energy 
eM is 

'M\ a495 /Qo Na742 



e« 0.0001855 {—) ^) . (53) 

If one multiplies this efficiency by the mass-energy M of the core, one gets that 
the total energy emitted in positrons is 



eM « 3.315 x W^Y ~ (9*™ 



M Q J \MJ 



ergs, (54) 



For stellar mass cores, this is 2-3 orders of magnitude smaller than the energy of 
gamma-ray bursts [Ruffini et al.(1998), Ruffini et al. (1999a), Ruffini et al. (1999b), 
Ruffini et al.(2000), Torres & Anchordoqui(2004)], essentially because the efficiency 
is 3-4 orders of magnitude less than unity. 
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Again, I should emphasize that all of these estimates and calculations give only 
very conservative upper limits on the efficiency and energy emitted, since they all 
assume that somehow one can hold the charge onto the surface of the collapsing core 
even when it is as large as the transition radius R t , which for £ — Qo/M > 10 -4 
would be significantly greater than the size of a neutron star. For such a core the 
electrostatic forces of repulsion on the excess protons on its surface would be more 
than 14 orders of magnitude greater than the gravitational attraction of the core for 
these protons, and it is hard to imagine that any forces sufficiently powerful to hold 
in the protons (such as nuclear forces at nuclear distances) could be effective in any 
star larger than a neutron star. And if one does take £o ~ 10~ 4 so that Rt ~ 50 km, 
somewhat larger than a neutron star size, then the upper limit on the efficiency of 
the pair production process would be only of the order of 2 x 10~ 7 , which is far too 
small to give a viable model for gamma ray bursts. 

6 Conclusions 

It does not seem to be possible to have astrophysical dyadospheres [electric fields 
larger than the critical value for Schwinger pair production, over macroscopic regions 
much larger in all directions than the regions of high fields that might conceivably 
be produced by individual heavy nuclear collisions, or much thicker than possible 
very thin high-field regions at the surface of a strange star [Usov(1998), Usov(2004), 
Usov, Harko, & Cheng(2005)] or neutron star [Ruffini(2006)]]. If, as is most plau- 
sible, charge carriers like protons are bound to an astrophysical object, such as a 
star or stellar core, primarily by gravitational forces, then the electric field cannot 
get within 13 orders of magnitude of the minimal dyadosphere values. (The excess 
charges will simply be ejected by the electrostatic repulsion when that exceeds the 
gravitational attraction. Then pair production rates from the macroscopic electric 
field, as opposed to that from collisions of individual particles, will be trillions of 
orders of magnitude below dyadosphere values and so will be completely negligible.) 

Even in what I consider to be the implausible scenario in which the excess charge 
carriers are bound by nuclear forces to a collapsing stellar core, I have shown here in 
a simple spherically symmetric model that the electric field has a very conservative 
maximum value that is more than a factor of 18 below the minimal dyadosphere 
value. Because the pair production rate is essentially exponential in the negative 
inverse of the electric field, the upper limit of the pair production rate in even 
this implausible scenario is more than 26 orders of magnitude below the minimal 
dyadosphere values. 

The idealized implausible scenario considered here, to give this very conservative 
upper bound on the pair production rate, had the maximal amount of charge some- 
how bound to the surface of an idealized stellar core with maximal initial charge 
that undergoes free fall collapse from radial infinity in the Schwarzschild metric 
(conservatively ignoring the fact that an actual astrophysical collapse would start at 
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finite radius and fall in slower, giving more time for discharge, and the fact that if 
the initial charge were maximal, the strong electric field would modify the geometry 
and also give electrostatic repulsion of the core surface, both of which would also 
slow down the collapse and lead to greater discharge and smaller electric fields). 
This scenario led to the maximal electric field (occurring when the freely collapsing 
core enters the event horizon of the black hole that would form) being less than 
5.5% of dyadosphere values. Although the pair production rate is always less than 
1CT 26 that of dyadosphere values, in this implausible scenario of having no other 
mechanism of discharging the core, there would be enough pair production to keep 
the electric fields always more than 18 times smaller than dyadosphere values. 

Although macroscopic astrophysical dyadospheres do not form in this example or 
in any other similar example that has been considered, in this idealized implausible 
scenario there is significant pair production (though at macroscopic astrophysical 
rates that are much, much lower than microscopic dyadosphere rates), and so I 
calculated what fraction of the total mass-energy M of the collapsing stellar core 
would be converted into pairs. Here I found that this efficiency, even under the highly 
idealized conditions of having maximal initial charge at such large radii that it seems 
inconceivable that the charge carriers could be sufficiently bound to such objects 
so much larger than neutron stars, is always much less than unity for collapsing 
objects with much less mass than three million solar masses: the efficiency is very 
conservatively bounded by 2x lCT 4 ^ M/M & . Therefore, even these idealized charged 
collapsing objects, unless they were enormously more massive than the sun, would 
not produce enough energy in outgoing charged particles to be consistent with the 
observed gamma ray bursts. 

It would of course be of interest to calculate the upper limits on the electric 
field and on the maximum pair production rate for models in which one relaxed the 
spherical symmetry. Although I would readily admit that I do not have a rigorous 
proof that the pair production rate cannot be higher then, the general arguments of 
Section 3 strongly suggest that it would be very surprising if it could be much larger. 
Therefore, the spherically symmetric model analyzed here lead me to conjecture that 
the maximal pair production rates achievable by astrophysical electric fields that are 
macroscopic in all directions are always less than roughly 10~ 26 that of hypothetical 
dyadosphere values. 

In conclusion, macroscopic dyadospheres almost certainly cannot form astrophys- 
ically, and the much weaker pair production rates that might occur, under highly 
idealized and implausible scenarios, do not seem sufficient for giving viable models 
of gamma ray bursts. 

For a shorter version of this work for a conference proceedings, see [Page(2006b)]. 
For more details, see [Page(2006a)]. 
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